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Abstract 


The paper is concerned with conservative solutions to the nonlinear wave equation 
utt—c{u){c{u)ux)^ = 0. For an open dense set of initial data, we prove that the solution 
is piecewise smooth in the t-x plane, while the gradient Ux can blow up along finitely many 
characteristic curves. The analysis is based on a variable transformation introduced in [7], 
which reduces the equation to a semilinear system with smooth coefficients, followed by 
an application of Thom’s transversality theorem. 


1 Introduction 


Consider the quasilinear second order wave equation 


utt - c{u){c{u)ux)x = 0 


t G [0,r], X G M. 


( 1 . 1 ) 


On the wave speed c we assume 


(A) The map c : M i—>• M_|_ is smooth and uniformly positive. The quotient c'{u)/c{u) is 
uniformly bounded. Moreover, the following generic condition is satisfied: 




( 1 . 2 ) 


Notice that, by (1.2), the derivative c'{u) vanishes only at isolated points. 
The analysis in [7, 3] shows that, for any initial data 


^(O,^) = uo{x), rii(0,x) = ui{x), 


(1.3) 


with uq G 77^ (M), ui G L^(R), the Cauchy problem admits a unique conservative solution 
u = u{t,x), Holder continuous in the t-x plane. We recall that conservative solutions satisfy 
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an additional conservation law for the energy, so that the total energy 

£{t) = ^ J[u^ + c^{u)ul]dx 

coincides with a constant for a.e. time t. A detailed construction of a global semigroup of 
these solutions, including more singular initial data, was carried out in [15]. 

In the present paper we study the structure of these solutions. Roughly speaking, we prove 
that, for generic smooth initial data (uo,rti), the solution is piecewise smooth. Its gradient Ux 
blows up along finitely many smooth curves in the t-x plane. Our main result is 

Theorem 1. Let the funetion u c{u) satisfy the assumptions (A) and let T > 0 be given. 
Then there exists an open dense set of initial data 

V C (c^{R) n X (^C2(M) nL2(R)^ 

sueh that, for {uo,ui) € V, the conservative solution u = u{t,x) of (1.1)-(1.3) is twiee con¬ 
tinuously differentiable in the complement of finitely many characteristie eurves 7 ^, within the 
domain [0,T] x M. 

For the scalar conservation law in one space dimension, a well known result by Schaeffer [17] 
shows that generic solutions are piecewise smooth, with finitely many shocks on any bounded 
domain in the t-x plane. A similar result was proved by Dafermos and Geng [ 8 ], for a special 
2x2 Temple class system of conservation laws. It remains an outstanding open problem to 
understand whether generic solutions to more general 2 x 2 systems (such as the p-system of 
isentropic gas dynamics) remain piecewise smooth, with finitely many shock curves. 

The proof in [17] relies on the Hopf-Lax representation formula, while the proof in [ 8 ] is based 
on the analysis of solutions along characteristics. In the present paper we take a quite different 
approach, based on the representation of solutions in terms of a semilinear system introduced 
in [7]. In essence, the analysis in [7] shows that, after a suitable change of variables, the 
quantities 


w = 2 arctan(ut + c(u)w 3 ;), z = 2 a,icta,n{ut — c{u)ux), 

satisfy a semilinear system of equations, w.r.t. new independent variables X, Y. See (2.16)- 
(2.20) in Section 2 for details. Since this system has smooth coefficients, starting with smooth 
initial data one obtains a globally defined smooth solution. To recover the singularities of the 
solution u of ( 1 . 1 ) in the original t-x plane, it now suffices to study the level sets 

{w{X,Y) = 7r}, {z{X,Y) = 7r}. (1.4) 

Since w and 2 ; are smooth, the generic structure of these level sets can be analyzed by tech¬ 
niques of singularity theory [ 2 , 9, 10, 14, 18], relying on Thom’s transversality theorem. One 
should be aware that, while the map {X,Y) i-)- {t,x,u,w, z) is smooth, the inverse map 
{t,x) I—> {X,Y) can have singularities. This variable transformation is indeed the source of 
singularities in the solution u = tt(t, x) of ( 1 . 1 ). 

The present work was motivated by a research program aimed at the construction of a distance 
which renders Lipschitz continuous the semigroup of conservative solutions of (1.1). Toward 
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this goal, one needs a dense set of piecewise smooth paths of solutions, whose weighted length 
can be controlled in time. In the final section of this paper we thus consider a 1-parameter 
family of initial data A {uq,u^), with A G [0,1]. We show that it can be uniformly 
approximated by a second path of initial data A i-)- (uq,u^), such that the corresponding 
solutions = u^{t,x) of (1.1) are piecewise smooth in the domain [0,T] x R, for all except 
at most finitely values of A G [0,1]. An application of this result to the construction of a 
Lipschitz metric will appear in the forthcoming paper [4]. 

The remainder of the paper is organized as follows. In Section 2 we review the variable change 
introduced in [7] and derive the semilinear system used in the construction of conservative 
solutions to (1.1). In Section 3 we construct families of smooth solutions to the semilinear 
system, depending on parameters. By a transversality argument, in Section 4 we show that for 
almost all of these solutions the level sets (1.4) satisfy a number of generic properties. After 
these preliminaries, the proof of Theorem 1 is completed in Section 5. Finally, in Section 6 
we prove a theorem on generic regularity for 1-parameter family of solutions. 

For the nonlinear equation (1.1), the formation of singularities in finite time was first studied 
in [11]. Based on the representations [7, 5], a detailed asymptotic description of structurally 
stable singularities is given in [6], for conservative as well as dissipative solutions. 

We conjecture that the regularity property stated in Theorem 1 should also hold for generic 
dissipative solutions of (1.1). However, in the dissipative case the corresponding semilinear 
system derived in [5] contains discontinuous terms, and smooth initial data do not yield 
globally smooth solutions. For this reason, the techniques used in this paper can no longer 
be applied. We remark that, at the present time, the uniqueness and continuous dependence 
of dissipative solutions to (1.1) has not yet been proved, for general initial data (uo,ui) G 
H^{R) X L2(R). 


2 Review of the main equations 


Consider the variables 


so that 


R = ut + c{u)ux , 
S = ut- c{u)ux , 


ut = 


R + S 


Ut = 


R-S 

2c 


For a smooth solution of (1.1), these variables satisfy 

Rt-cR^ = i{R^-S^), 

St + cS, = ^^{S^-R^). 

In addition, R^ and S'^ satisfy the balance laws 

{R%-{cR^)t = f^{R^S-RS^), 

{S^)t + {cS% = f^{S^R-SR^). 

As a consequence, for smooth solutions the following quantity is conserved: 

E = = - - -. 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 

(2.4) 

(2.5) 
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One can think of B? and as the energy densities of backward and forward moving waves, 
respectively. Notice that these are not separately conserved. Indeed, by (2.4) energy can be 
exchanged between forward and backward waves. 

It is well known that, even for smooth initial data, the quantities ut,Ux can blow up in finite 
time [11]. To deal with possibly unbounded values of i?, 5, following [7] it is convenient to 
introduce a new set of dependent variables: 

w = 2 arctani?, 2 : = 2arctan5. ( 2 - 6 ) 

Using (2.3), we obtain the equations 


2 , 

d - 52 


Wt CWx = 

2 c l + i?2 ’ 

(2.7) 

2 , 

c' 


Zt + CZx = g2^St +cSx) 

2c 1 + 52 ■ 

(2.8) 



Figure 1: Characteristic curves. As new coordinates of the point {t,x) we choose the values (AT, T) = 
(a:“ (0, t, x), — a;“''(0,t, x)). 

We now perform a further change of independent variables (Fig. 1). Consider the equations 
for the backward and forward characteristics: 

x~ = — c{u ), x'^ = c{u) , (2-9) 

where the upper dot denotes a derivative w.r.t. time. The characteristics passing through the 
point (t, x) will be denoted by 

s !-)■ x~{s,t,x), s x~^{s,t,x), 

respectively. As coordinates {X, Y) of a point {t, x) we shall use the intersections of these 
characteristics with the x-axis, namely 

X = x~{0,t,x), Y = —x~^{0,t,x). (2-10) 

Of course this implies 

Xt-c{u)Xx = 0 , Yt + c{u)Yx = 0 , ( 2 - 11 ) 
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— 0, + (c^)x — 0. (2-12) 

For any smooth function /, using (2.11) one finds 

ft + cU = fxXt + fYYt + cfxX^ + cfyY, = {Xt + cX,)fx = 2cX,fx, 


ft-cU = fxXt + fyYt-cfxX^-cfYY, = iYt-cY,)fY = -2cY,fY. 


We now introduce the further variables 

l + i ?2 


P = 




Notice that the above definitions imply 

1 p o w 


l + i ?2 = 2 


q = 


-1 


1 + 52 
-w 


rr:^ = 2 


(2.13) 

(2.14) 

(2.15) 


Starting with the nonlinear equation (1.1), using X,Y as independent variables one obtains a 
semilinear hyperbolic system with smooth coefficients for the variables u,w, z,p,q, namely 


J 

j ux 

sin 117 ^ 

~ 4c F) 


1 

[ uy 

- ^0 

4c q ’ 


Wy 

c’ 

8? 

T (cos Z — COS W 

)q 

Zx 

_ c' 

8c? 

T (cos W — COS Z 

)p 

py 

c' 

8i? 

(sin z — sin w) 

pq 

qx 

c' 

8i? 

(sin w — sin z) 

pq 


(2.16) 


(2.17) 


(2.18) 


The map {X,Y) i-). {t,x) can be constructed as follows. Setting f = x, then / = t in the two 
equations at (2.13), we hnd 


r c = 2cX^ XX , 

I -c = -2cY^xy, 

respectively. Therefore, using (2.15) we obtain 


1 = 2cXx tx , 
1 = - 2 cya ; ty , 


XX = 

1 

2Xj; 

(1+cos w) p 
~ 4 

Xy = 

1 

(l+cos z) q 
~ 4 

f ix = 

1 

2cX^ 

_ (H-COS117)p 

~ 4c ’ 

i W = 

1 

(l+cos z) q 

-ScVc 

~ 4c 


(2.19) 


( 2 . 20 ) 


Given the initial data (1.3), the corresponding boundary data for (2.17)-(2.20) can be deter¬ 
mined as follows. In the X-Y plane, consider the line 


70 = {x + y = o} c 
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parameterized as x i-)- {X{x), 5^(x)) = (x, —x). Along 70 we can assign the boundary data 
(u, W, z, p, q) by setting 


_ _ . . ( w = 2 arctan i?( 0 , x), ( p = l + i?^( 0 , x), 

u uo[x), I ^ _ 2 arctan ^(O, x), \q = 1 + S'^(0,x), 

at each point (x, —x) G 70 . We recall that, at time t = 0, by (1.3) one has 

i?( 0 ,x) = {ut + c{u)ux){0,x) = ui{x) + c{uo{x))uo,xix), 

S{0,x) = (ut - c{u)ux)i0,x) = ui{x) - c{uo{x))uo,xix). 


( 2 . 21 ) 


Remark 1. Since the semilinear system (2.17)-(2.20) has smooth coefficients, for smooth 
initial data all components of the solution remain smooth on the entire X-Y plane. As proved 
in [7], the quadratic terms in (2.18) (containing the product pq) account for transversal wave 
interactions and do not produce finite time blow up of the variables p, q. Moreover, if the values 
of p,q are uniformly positive along a fine {X + Y = k}, then they remain uniformly positive 
on compact sets of the X-Y plane. Throughout this paper, we always consider solutions of 
(2.17)-(2.20) where p,q > 0. 

By expressing the solution u{X,Y) in terms of the original variables (t,x), one obtains a 
solution of the Cauchy problem (1.1)-(1.3). Indeed, the following was proved in [7]. 

Lemma 1. Let {u,w, z,p,q,x,t) be a smooth solution to the system (2.16)-(2.20), with 
p,q> 0. Then the set of points 

{{t{X,Y),x{X,Y),u{X,Y)); (A,y)€M 2 | (2.22) 

is the graph of a conservative solution to the variational wave equation (1.1). 

We observe that, while the functions 

{X,Y) u{x,Y), (A,y) ^ A(A,y) = (t(A,y), x(A,y)) (2.23) 

are globally smooth, the map A : may not have a smooth inverse. Indeed, A may not 

even be one-to-one. Therefore, the solution u{t,x) = u{A~^ {t, x)) can fail to be smooth. This 
happens precisely at points where the Jacobian matrix DA is not invertible. By (2.19)-(2.20), 
singularities occur when cosw = — 1 or cosz = — 1 . 


Remark 2. The system (2.17)-(2.20) is invariant under translation by 27r in w and z. We 
can thus think of w,z as points in the quotient manifold T = R/27rZ. Throughout the 
following we take advantage of this fact and regard a solution of (2.17)~(2.20) as a map 
{X, Y) i-> («, w, z,p, q, X, t) from R^ into RxTxTxRxRxR. Observe that we have the 
implications 

cos w > — 1 , 


W ^ TT 
Z TT 
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cos z > — 1 . 


(2.24) 



Remark 3. In general, many distinct solutions to the system (2.16)-(2.20) can yield the same 
solution u = u{t,x) of (1.1). 

Indeed, let {u,w, z,p,q,x,t){X,Y) be one particular solution. Let : M I— )■ R be two 
bijections, with cp' > 0 and > 0. Introduce the new independent and dependent variables 
{X,Y) and {u,w,z,p,q,x,t) by setting 


A = (/)(A), 

Y = i;{Y), 

(2.25) 

u, w, z, 

x,i){X,Y) 

= {u,w,z,p,q,x,t){X,Y), 

(2.26) 

1 

r p{x,Y) 

= p{X,Y)-cj)'{X), 



1 

1 ^ ^ 


(2.27) 

1 

1 QiX,Y) 

= q{X,Y)-i;'{Y). 



Then, as functions of {X,Y), the variables {u,w, z,p, q, x,i) provide another solution of the 
same system (2.16)“(2.20). Moreover, by (2.26) the set 

[{i{X,Y),x{X,Y),u{X,Y)y, (X,y)GR2} (2.28) 

coincides with the set (2.22). Hence it is the graph of the same solution u of (1.1). One 
can regard the variable transformation (2.25) simply as a relabeling of forward and backward 
characteristics, in the solution u. A detailed analysis of relabeling symmetries, in connection 
with the Camassa-Hom equation, can be found in [13]. 

For future reference we observe that 


w^{X,Y) = wx{X,Y)-cl>'{X), 

Wj^^{X,Y) = wxx{X,Y) ■ [<j,'{X)]^ + wx{X,Y) ■ ^"(X). 

In particular, one has the equivalences 

w^iX,Y) = 0 

ZyiX,Y) = 0 

w^^){X,Y) = ( 0 , 0 ) 

{Zy, Zyy){X,Y) = (0,0) 

2.1 Compatible boundary data 

More generally, instead of (2.21) we can assign boundary data for the system (2.13)“-(2.18) on 
a line 7 = {X + Y = k}. Namely: 

P{s ), 
q{s) , 
(2.30) 


u(s, K — s) = u{s) , 


w{s, K — s) 
z{s, K — s) 


= w[s) 
= z{s) 


P{S,K- S) 

q{s,K- s) 


wxiX,Y) = 0, 
zy{X,Y) = 0, 
{wx,wxx){X,Y) = (0,0), 
{zy,zyy){X,Y) = (0,0). 


(2.29) 
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(2.31) 


for suitable smooth functions u,w, z,p, q. If both identities in (2.16) hold, then 


d . d , . sin w 

—u{s) = —u{s,K-s) = ux-UY = — — p- 
as as 4c 

The boundary data should thus satisfy the compatibility condition 


sinz 

4c 


q- 



sinr(;(s) 

4c(u(s)) 


p{s) 


sin z{s) 
4c(m(s)) 


q{s) 


(2.32) 


As remarked earlier, the system (2.16)-(2.18) is overdetermined. Indeed, the function u = 
u{X, Y) could be recovered by either one of the identities in (2.16). We now prove that, if the 
compatibility condition (2.32) holds, then any smooth solution satisfying one of the identities 
in (2.16) satisfies the other as well. 


Lemma 2. Let u,w, z,p,q he smooth functions on which satisfy (2.17)-(2.18) together with 
the boundary conditions (2.30) along the line 7 = {X + Y = k}. Assume that the compatibility 


condition (2.32) holds. 

Then 

one has 




uy - 

sinz 

Ac{u) ^ 

for all (A, Y) € 

(2.33) 

if and only if 

ux -- 

siiiw 

4c(ci) ^ 

for all (A,y)GR2. 

(2.34) 


Proof. Consider the smooth, strictly increasing function 

pU 

<l>(n) = / 4:c{s)ds. 

Jo 

Observe that the identities (2.33), (2.34) are equivalent respectively to 

^{u)y = sinz-( 7 , ^iu)x = sinro-p. (2.35) 


Assume that (2.33) holds. Then 

d>(n(A,y)) = ^{u{X, K- X)) + [ [sinz-g](A,s)(is. (2.36) 

Jk-X 

Differentiating w.r.t. X, and using the first equations in (2.17)-(2.18) together with the com¬ 
patibility condition (2.32) we obtain 

^{u)x{X,Y) = ^'{u)-[ux-UY]iX,K-X) +[sin z-q]{X,K-X) 


+ / [cos z ■ zxq + sinz ■ qx]{X, s) ds 
Jk-Y 


= [sinrc • p](A, K — A)-|- 


" c'(u) 


_Y L8c2(n) 


(1 -|- sin(z -|- w))pq (A, s) ds 


fV 


Ik-Y 


d 


= [since • p](A, K — A)-|- j ^^[sinu; • p](A, s) ds = [since • p](A, A). 

(2.37) 



We have thus proved that the second identity in (2.35) holds. This is equivalent to (2.34). 
The converse implication is proved in the same way. □ 


Next, consider initial data for t, x, on the curve 7 = {X + Y = k}, say 

x{s,K — s) = x{s), t{s,K — s) = t{s). (2.38) 

Using (2.19)-(2.20) we derive the compatibility conditions 



(1 + cos q{s))p{.s) + (1 + cos z{s))q{s) 

(2.39) 

4 

d , 

(1 + cosw{s))p{.s) — (1 + cos z(s))g(s) 

(2.40) 

4c(u(s)) 


Lemma 3. Let {u,w, z,p,q){X,Y) be a solution of the system (2.16)-(2.18). Then there 
exists a solution {t,x){X,Y) of (2.19)-(2.20) with boundary data (2.38) if and only if the 
eompatibility conditions (2.39)-(2.40) are satisfied. 


Proof. 1. Assume that the equations (2.19)-(2.20) are satisfied for all {X,Y) G M^. In 
particular, they are satisfied along the curve 7 = {X + Y = k}. This implies 


— x(s) = —x(s,K — s) = [xx — a:v](s, K — s) = 
ds ds 


(1 + cos w)p + (1 + cos z)q 


where the right hand side is evaluated at {X,Y) = (x,k — s). Hence (2.39) holds. The identity 
(2.40) is derived in the same way. 


2. Next, assume that the compatibility conditions (2.39)-(2.40) are satisfied. To prove that 
(2.19) admits a solution, it then suffices to check that the differential form 

(1 + cos w)p _ (1 + cosz)g 


is closed. This is true because 
(1 + cos w)p 


1 + cos w 


J y 


sintc 

-^WYP+ ^ 


■PY 


sm w c' , , 1 + cos w c' 


4 z-cos w)pq+ ^ 


8c2 


32c2 


p n 


(1 + COS w) sin z — (1 + cos z) sin w 

pq = 


(sin z — sin w) pq 
(1 + cos z)q 


(2.41) 

Similarly, to prove that (2.20) admits a solution, it suffices to check that the differential form 

(1 + cos w)p (1 + cosz)g 


4c 


4c 
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is closed. This is true because 


(1 + cos w)p 
4c 


sintc 


4c 


-wy p 


1 + cos w , 


4c2 


cuy p + 


1 + cos w 
4c 


■PY 


sin wc 1 + cos w , sin z 1 + cos w c 

^(cosz - cosw)pq --c -^m+ -- 


4c 8c2 


32c3 


(1 + cos w) sin z + {1 + cos z) sin w 


4c 4c 

(1 + cos z)q 


pq = 


4c 




sin w) pq 


(2.42) 

□ 


Remark 4. Let a solution (u, w, z,p, q) of (2.16)-(2.18) be given. If we assign the values of t, x 
at a single point {Xq, Yq)) then by the compatibility conditions (2.39)-(2.40) and the equations 
(2.19)-(2.20) the functions t{X,Y), x{X,Y) are uniquely determined for all {X,Y) G R^. 
Choosing different values of t, x at the point (Xq, Yq) we obtain the same solution u = u{t, x) 
of ( 1 . 1 ), up to a shift of coordinates in the t,x plane. 


3 Families of perturbed solutions 

Let a point (Xq, Yq) be given and consider the line 

7 = {(X,y); X + Y = k], k = Xq + Yq. (3.1) 

We can then arbitrarily assign the values of w, z,p, q at every point (X, Y) G 7 . Moreover, we 
can arbitrarily choose the values of u,x,t at the single point (Xo,Xo)- In turn, these choices 
uniquely determine functions u, x, t on 7 which satisfy the compatibility conditions (2.32) and 
(2.39)-(2.40). 

Based on this observation, we can construct several families of perturbations of a given solution 
of (2.16)-(2.18). The main goal of this section is to prove 


Lemma 4. Let the assumption (A) hold. Let {u,w, z,p, q) be a smooth solution of the system 
(2.16)-(2.18) and let a point {Xq,Yo) G R^ be given. 

(1) If {w,wx,wxx){Xq,Yo) = (tTjO,0), then there exists a 3-parameter family of smooth 
solutions {u^,w^, z^,p^,q^) of (2.16)-(2.18), depending smoothly on 9 G R^, such that 
the following holds. 

(i) When 0 = 0 G R^ one recovers the original solution, namely {u^,w^,z^,p^,q^) = 
{u,w,z,p, q). 

(ii) At the point (Xq, Yq), when 6 = 0 one has 

rank , Wxx) — 3. (3-2) 

(2) If {w, z,wx){Xq,Yq) = (tt, 7 r,0), then there exists a 3-parameter family of smooth solu¬ 
tions {u^,w^, z^,p^,q^) satisfying (i)-(ii) as above, with (3.2) replaced by 

rank He(rc® , z^ , Wx) = 3. (3.3) 
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(3) If (w,wx,c'(u))(Xq,Yq) = (tt, 0,0), then there exists a 3-parameter family of smooth 
solutions ,w^, ,p^,q^) satisfying (i)-(ii) as above, with (3.2) replaced by 


iBvak D0{w^ , Wxi c'{u^)) = 3. 


(3.4) 


For example (3.2) means that we can construct perturbed solutions, depending on parameters 
01,02,03, such that the Jacobian matrix 


/ 


D0{w^, wx, '‘^xx) — 








■W^wxx 

m^wxx 

m^wxx 


\ 


computed at 0 = 0, has full rank at the point {Xq, Yq)- 

3.1 Proof of Lemma 4 


(3.5) 


Let {u,w, z,p,q) be a solution of the system (2.16)-(2.18). Given the point (Xo,lo)) 
consider the line 7 in (3.1) and let {u,W,z,p,q) be the values of the solution along 7 , as in 
(2.30). 

For future use, we compute the values of wx,wxx at the point (Xo,lo)- At any point 
(s, /t — s) € 7 we have 


wx — VJy = w'(s) 


Zx - ZY = ^ (s) 


qx-qr = q'{s). 


Here and in the sequel, a prime denotes derivative w.r.t. the parameter s along the curve 7 . 
Using (2.17)-(2.18) we obtain 


wx{Xo,Yq) = w' + ^ (cos z - cos w)q, 

Sc^{u) 

d (w) 

Zy{Xo,Yo) = - ^ ^ (cosIZJ - cos z)p , 

Sc^[u) 

d (w) 

gy(Xo,Fo) = -9' +x-774(sinIZJ-sinz)^, 

8c^{u) 

where all terms on the right hand sides are evaluated at s = Xq. 

A further differentiation yields 
(f _ d 

^^^^(s) = -^J(wx{s,n-s)- wy{s,k- s)] = [wxx + wyy -‘Iwxy]{s,k - s). 


(3.6) 

(3.7) 

(3.8) 
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Using (2.16)-(2.18) together with (3.6)~(3.8) we obtain 
yJYX i^o, ^o) 

= ux{cosz-cosw)q + -^{wxsmW-zx sin z)q + -^ {cos z-cosW)qx 

= (^)'^(cosz-cos«J)m 

+ ^ + ^(cosz — cosw)q) sinu; — -^{cosw — cosz)psinz| q 

+ (^)^ (cosz — coswJ) (sinifj — sin z)pq 


= fi 


and 


(3.9) 


wyy{Xo, Uq) 

= ^ sJ \^) ) (cosis — cos uJ) g ^ (iZJy sinla — zy sinz) g + ^ (cosz — cosIZJ) gy 

= (^)'^(cosz-cosIZJ)g2 

+ ^ |^(cosz — cosW)qsinw — { — z' + -^{cosw — cos'z)p) sinz| q 
+ ^ (cos z — cosw) { — q' + ^(sinlZJ — sin z)pq) 


= /2. 


Hence 

wxx{Xo, Uo) = w" + 2/i — /2 . 


(3.10) 

(3.11) 


We now construct families {u^, q^, z^,p^, q^) of perturbations of the data (2.30) along the curve 


7 , so that at the point (Xo,hb) the matrices in (3.2)~(3.4) have full rank. These perturbations 

will have the form 



( w\s) = w{s) + YPl=iOiWi{s), 

( /(s) = p{s) + J2Li^iPiis), 

(3.12) 

\ z\s) = %s) + Y.l=iOiZi{s), 

\ Q^{s) = q{s) + Td=i9iQiis), 


for suitable functions lUj, Zj, Pj, Qj G . 

Moreover, at the point s = Xq we set 


u^iXo) = 

u(Xo)+ e^Ui. 

(3.13) 


i=l,2,3 


In turn, the above definitions together with the compatibility conditions (2.32) determine the 
values of viP (s) for all s € M. In particular, for each 0 E we obtain a unique solution of the 
semilinear system (2.16)~(2.18). 
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We observe that the functions Wi, Zi, Pi,Qi can be chosen arbitrarily. Hence at the point 
s = Xq, 9 = 0, we can arbitrarily assign all derivatives 

d ^ d d^ _g d d^ _g d d^ 

dO ds^^ ’ dO ds^^ ’ dO ds^^ ’ dO ds^^ ’ 
with /c = 0,1, 2,... Moreover, we can arbitrarily choose the quantity ■^u^{Xo), while all higher 
order derivatives with k > 1, are then determined by the compatibility condition 

(2.32). 


1. To achieve (3.2), we choose perturbations i = 

Jacobian matrix of first order derivatives w.r.t. 6 *i, 6 * 2 , 6 * 3 , computed at s 
given by 


Da 


Da 


/ u \ 


( ° 

0 


w 


1 

0 

0 

1 


0 

0 

0 

z' 


0 

0 

0 

w' 

= 

0 

1 

0 

w" 


0 

0 

1 

p 


0 

0 

0 

q 


0 

0 

0 

\ 1' 


V 0 

0 

0 / 

i (3.11) 

this yields 



w 

\ 

/I 

0 

0 

wx 

= 

* 

1 

0 

Wxx 

/ 

\ * 

* 

1 


1,2,3, so that the 
= Xq and 6 * = 0 , is 


Notice that, for the third family of perturbations (corresponding to the third column), the 
first order variations of fi and /2 in (3.9)-(3.10) both vanish at (Jfo,lo). This achieves (3.2). 


2. To achieve (3.3), we choose perturbations 


s = Xq and 0 = 0 one has 

( u \ 


f 0 

0 

0 \ 


w 


1 

0 

0 

De 

'z 


0 

1 

0 

w' 

— 

0 

0 

1 


p 


0 

0 

0 


\ Q / 


V 0 

0 

0 / 

At the point (Ao,Fb)) by (3.6) this yields 





/ 

^ w ^ 

\ 

/ 1 

0 

0 \ 

De 

z 

= 

0 

1 

0 

\ 

V wx J 

/ 

\ * 

* 

1/ 


Hence (3.3) holds. 


1,2,3, so that at 


(3.14) 


(3.15) 


3. Finally, we construct three families of perturbations satisfying (3.4). If at (Xo,Fb) we have 
c'(u(Xo,To)) = 0, then the assumption (A) implies 

c"(n(Ao,yo)) / 0. (3.16) 


13 



To achieve (3.4), we choose three families of perturbations such that, at s = Xq and 9 = 0, 


/ 

u 

\ 


(° 

0 

1 \ 


w 



1 

0 

0 


'z 



0 

0 

0 


w' 



0 

1 

0 


p 



0 

0 

0 

V 

1 

) 


V 0 

0 

0 / 


At the point (Ao,lo)) by (3.6) and the first equation in (2.17), this yields 


De wx 
\ c'{u) 



0 

1 


0 

* 


'{u) 


This achieves (3.4). 


(3.17) 


(3.18) 

□ 


4 Generic solutions of the semilinear system 


In this section we study smooth solutions to the semilinear system (2.17)~(2.20), determining 
the generic structure of the level sets {{X,Y ); w{X,Y) = vr} and {{X,Y); z{X,Y) = n}. 


Lemma 5. Let the function u i-)- c{u) satisfy the assumptions (A) and consider a compact 
domain of the form 

r = {{X,Y); |A| + |y| < m}. (4.1) 

Calls the family of all solutions to the system (2.16)-(2.18), withp,q > 0 for all {X,Y) G 
Moreover, call S' C S the subfamily of all solutions {u, w, z,p, q) such that, for {X, Y) G T, 
none of the following values is attained: 


{w,wx,wxx) 

II 

o 

o 

(4.2) 

{z,zy,zyy) 

II 


II 

T 

(7r,7r,0), 

(4.3) 

II 

o 

T 

j) 

(7r,7r,0), 


{w,wx,c'{u)) 

II 

o 

(4.4) 

{z,ZY,c'{u)) 

II 

o 



Then S' is a relatively open and dense subset of S, in the topology induced by C^(r). 

Some words of explanation are in order (Fig. 2). Asking that the values in (4.2) are never 
attained is equivalent to the implications 


w = TT and Wx = 0 
z = IT and zy = 0 


wxx 7 ^ 0, 

Zyy 7^ 0. 


(4.5) 
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Figure 2: Two level sets {w = tt} and {z = tt}, in a generic solution of (2.16)-(2.18). At Pi, P 2 one 
has w = TT, wx = 0 while the generic conditions imply wy ^ 0, wxx 7 ^ 0. At the points Qi, Q 2 where 
the two singular curves cross, by (2.17) one has wy = zx = 0, while the generic conditions imply 
Wx 7 ^ 0, zy 7 ^ 0. Hence the two curves have a perpendicular intersection. 

Writing the level curves in the form {w{X,Y) = vr} = {Y = ip{X)'\ and {z(X,Y) = tt} = 
{X = ^(y)}, this imposes some restrictions at the points where = 0 or = 0. 

Asking that the values in (4.3) are never attained is equivalent to the implication 

[w = TT and z = tt] => [wx 7 ^ 0 and zy 7 ^ 0]. (4.6) 

This imposes restrictions at points where two level curves {w = tt} and {z = tt} cross each 
other. 

Finally, the lemma states the existence of a perturbed solution such that values (4.4) are never 
attained. To understand the meaning of this condition, consider a solution which never attains 
any of the values in (4.3)-(4.4). In this case, by (2.24) the conditions w = tt and wx = 0 
together imply 

(w) 

= o ^ ^ 0 - 

8c^[u) 

This is equivalent to the implication 

W = TT {wx,Wy) 74 (0,0). 

By the implicit function theorem, the level set {w = tt} is then the union of regular curves in 
the X-Y plane (restricted to the domain F). Similarly, the level set {z = tt} will be a union 
of regular curves. 

We shall give a proof of Lemma 5, using Lemma 4 together with Thom’s transversality theo¬ 
rem. For readers’ convenience, we first review some basic definitions [2, 12, 18]. 


Definition (map transverse to a submanifold). Let / : X 1—7 T be a smooth map of 
manifolds and let IF be a submanifold of Y. We say that / is transverse to VF at a point 
p € X, and write / fflplF, if 
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Figure 3: A generic solution u = u{t,x) of (1.1) with smooth initial data remains smooth outside 
finitely many singular points and finitely many singular curves where Ux —> ±oo. The curves where u 
is singular are the images of the curves where re = tt or z = tt in Fig. 2, under the map (t, x) = A(X, Y) 
at (2.23). Here pi = A(Pi) are points where singular curves originate or terminate, while qj = ^{Qj) 
are points where two singular curves cross. 


• either f{p) ^ W, 

. or else f{p) G W and Tf(p)Y = (df)p(TpX) + 

Here TpX denotes the tangent space to X at the point p € X, while TgY and TgW denote 
respectively the tangent spaces to Y and to W at the point q £ W C Y. Finally, {df)p : 
TpX i-G Tj(^p^Y denotes the differential of the map / at the point p. 

We say that / is transverse to W, and write / ffl W, if / fflplF for every p £ X. 

In the special case where W = {y} consists of a single point, / fflIF if and only if y is a regular 
value of /, in the following sense. 


Definition (regnlar value). Let / : X i—F be a smooth map of manifolds. A point y £Y 
is a regular value if, for every p £ X such that f{p) = y, one has 

Tf^p)Y = {df)p{TpX). 


Transversality Theorem. Let X, 0, and Y he smooth manifolds, W a submanifold ofY. 
Let 9 fP he a smooth map which to each 6 £ Q associates a function fP £ C°°{X,Y), and 
define <I> : X x 0 i-> F by setting ^{x,6) = (fP{x). // rjl IF then the set {6 £ Q,; (fP t\\W} 
is dense in 0. 

For a proof, see [2, 12], 

4.1 Proof of Lemma 5. 

1. We shall use the representation 

S' = SinS 2 ns^nSins^r\SQ, (4.7) 
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where Si,... ,Sq C <S are the families of solutions for which one of the six values listed in 
(4.2)-(4.4) is never attained on F. For example, 5i is the set of all solutions such that 


{w,wx,wxx){X,Y) / (7r,0,0) 

for all (A, Y) G F, 

(4.8) 

while Sq is the set of all solutions such that 

{z,zy,c'{u))(X,Y) 7 ^ (7r,0,0) 

for all (A, Y) G F. 

(4.9) 


Since F is a compact domain, it is clear that each Si is a relatively open subset of S, in the 
topology of C^(F). In the remainder of the proof we will show that each Si is dense on S. 

2. Let {u,w, z,p,q) be any solution of (2.16)-(2.18), with p,q > 0. By a smooth approx¬ 

imation of the data along the line 70 = {X + Y = 0}, it is not restrictive to assume that 
u,w,z,p,q G We begin by looking at the first condition in (4.2). 

Given any point {Xq, Yq) G F, two cases can occur. 

CASE 1: {w,wx,vjxx){Xo,Yq) ^ (tt, 0,0). In this case, by continuity, there exists a neigh¬ 
borhood M of (Ao, Yq) in the X-Y plane where we still have {w, wx,wxx} 7 ^ (tt, 0, 0). 

CASE 2: {w,wx,wxx){Xq,Yq) = ( 7 r, 0 ,0). In this case, by Lemma 4 we can find a 3- 

parameter family of solutions {u^,w^ , z^,p^,q^) such that the 3x3 Jacobian matrix of the 
map 

(01,02,93) ^ (w^{X,Y), w<^x{X,Y),w$,x{X,Y)'^ (4.10) 

has rank 3 at the point (Xq, Yq), when 9 = 0. By continuity, this matrix still has rank 3 on a 
neighborhood M of (Ao,Fo)) for 9 sufficiently close to zero. 

We now choose finitely many points (Xi,Yi), i = 1,... ,n, such that the corresponding open 
neighborhoods Mf^XiXi) cover the compact set F. Call nj the cardinality of the set of indices 

X = {i; {w,wx,wxx)(Xi,Yi) = (7r,0,0)} (4.11) 

so that CASE 2 applies, and set N = 3ni. 

3. Let U D F be an open set contained in the union of the neighborhoods M(^Xi,Yi): ^'^icl call 

Bi; = {9 € ; |6*| < e} the open ball of radius e in . 

We shall construct a family (u^,w^ , z^,p^,q^) of smooth solutions to (2.16)"(2.18), such that 
the map 

iX,Y,9) ^ (w\X,Y), w<^x{X,Y), w<^xx(X,Y)') (4.12) 

from Q X into has (tt, 0, 0) as a regular value. Toward this goal, we need to combine 
perturbations based at possibly different points (Xi,Yi) into a single A-parameter family of 
perturbed solutions. 

Let (u, w, z,p, q)(X,Y) be a solution to the system (2.16)-(2.18). For each A: = 1,..., A, let a 
point (Afc, Yk) be given, together with a number G M and functions W^, Z^, P^, G C^(M). 
By the previous analysis, a 1-parameter family of perturbed solutions to (2.16)~(2.18) is then 
determined as follows. For |e| < sufficiently small, let 

(^i^^c^ 2 ;^/,g^) = ^l(u,w,z,p,q) (4.13) 
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be the unique solution of (2.16)-(2.18) with data assigned on the line 7 ^ = {X + Y = + 

by setting 

u%Xk,Yk) = u{Xk,Yk) + eUk , 

while for {X,Y) G 7 ^ 

= w + sWk , = z + eZk , = p + sPk , = q + eQk ■ 

Given ( 6 * 1 ,, On), a perturbation of the original solution {u, w, z,p, q) is dehned as the com¬ 
position of N perturbations: 

{v!^,w\z\p^,q^) = o • • • o (4.14) 


4. At each point (Xi,Yi) with i G I, we can apply Lemma 4 and obtain three 1-parameter 
families of perturbed solutions so that the Jacobian matrix (4.10) has rank 3 on Af(Xi,Yi), for 
all 6 small enough. 

Combining all these perturbations, we obtain an iV-parameter family of solutions such that 
the value (tt, 0,0) is a regular value for the map (4.12), from Q x into M^. 

By the transversality theorem, for a.e. 9 the value (tt, 0,0) is a regular value for the map 
{X,Y) {X,Y), w^{X,Y), w^j^{X,Y)'j from Q, into M^. Since Q, has dimension 2, 

for a.e. 9 the corresponding solution ,w^, z^,p^,q^) has the property that 

{w^{X,Y), w^xiX,Y), w^xx(X,y)) ^ (vr,0,0) 

for all (X, Y) G F. This proves that the set <Si of solutions for which (4.8) holds is dense on S. 

5. Repeating the above construction, we obtain that each i = 1,..., 6 , is a relatively open, 
dense subset of S. By (4.7), the intersection S' is is a relatively open, dense subset of S. □ 


5 Proof of Theorem 1. 

Consider the product space 

U= (c^{R)nH\R)'^ X (c‘^{R)nL‘^{R)^ (5.1) 

with norm 

11(^0,= ll^ollcs + ll^^olUi + ||lil||c 2 + II«i||l 2 • 

Given initial data (ho, hi) G U, consider the open ball 

B 5 = |(uo,ui)gW; ||(uo,ui) - (ho,hi)||^ < jj. (5.2) 

Theorem 1 will be proved by showing that, for any (ho, hi) G U there exists a radius J > 0 and 
an open dense subset V Bs, with the following property: For every initial data (no, ni) € V, 
the conservative solution n = u{t,x) of (1.1)-(1.3) is twice continuously differentiable in the 
complement of finitely many characteristic curves 7 ^, within the domain [0,T] x M. 
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1. Let (tio,wi) ^ U he given. By the definition of the space U in (5.1), as |x| —>■ oo we have 

uo{x) 0, uq^x{x) —>■ 0, ui{x) —)■ 0. (5-3) 

Hence the corresponding functions R,S in (2.1) satsfy 

i?(0, x) —7> 0, 5(0, x) ^ 0. 

From (2.3), it follows that the functions ii, 5 remain uniformly bounded on a domain of the 
form {(t,x); t G [0,T], |x| > r}, for r sufficiently large. More generally, we can choose 

5 > 0 such that, for every initial data {uq, ui) G Bs, the corresponding solution u{t, x) remains 

twice continuously differentiable on the outer domain 

|(t,x); tG[0,r], |x|>/)|, (5.4) 

for some p > 0 sufficiently large. Its singularities can thus occur only on the compact domain 

[0,r] X [-P, p\. 

The subset D d is now dehned as follows. (no,ui) G H if {uo,ui) G B^ and moreover, 
for the corresponding solution {u,w, z,p,q) of (2.16)-(2.20) with boundary data (2.21), the 
values (4.2)-(4.4) are never attained, for any {X,Y) such that 

{t{X,Y),x{X,Y)) G [0,T] X [-P, p]. (5.5) 

It is important to observe that, by (2.29), the above condition is independent of the relabeling 
(2.25). 

2. For any {uo,ui) G Bs we now consider the corresponding solution {t,x,u,w, z,p,q) of the 
system (2.16)”(2.20), with boundary data as in (2.21). Let A be the map at (2.23) and let F 
be the square with side 2M in the X-Y plane, as in (4.1). 

By choosing M large enough, and by possibly shrinking the radius 6, we can achieve the 
inclusion 

[0,T]x[-p,p] C A(r), (5.6) 

for every (uo,ui) G Bs- 

3. We begin by proving that V is open, in the topology of x C^. Indeed, consider initial 
data (uo,ui) G V and let (uq, u(')i,>i be a sequence of initial data converging to (uo,ui). 

Assume, by contradiction, that (uq,u^) ^ V for all iz > 1. To fix the ideas, let (X'^,Y‘') be 
points at which the corresponding solutions {u'^ ,w^, z^,p^,q'^) satisfy 

{w'^,w\,w\x){X\Y-) = (7r,0,0), (r,x")(X^y'^) G [0,T] x [-p, p], (5.7) 

for all V >1. By (5.6), since the domain F in (4.1) is compact, by possibly taking a subsequence 
we can assume {X'^,Y'') — {X,Y). By continuity, this implies 

{w,wx,wxx)iX,Y) = (tt, 0, 0), {t,x){X,Y) G [0,T] x [-p, p], 

contradicting the assumption (uo,ui) G V. 
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The other cases in (4.2)~(4.4) are handled in the same way. This proves that P is open. 


4. Next, we claim that D is dense in B^. Indeed, let {uo,ui) € Bs be given. By an arbitrarily 
small perturbation (measured in the norm of Z4), we can assume that uq,ui € C°°. 

Using Lemma 5, we can construct a sequence of solutions {u^,w'^, z'^,p''', , x'^,t^) of (2.16)- 

(2.20) such that: 


(i) For every bounded set 14 C and any k >1, the norm of the difference satisfies 


lim 

I/—>-oo 


w^ — Wj — Zj — Pj — Qj — X, —i)|| ^ = 0* (5-8) 


(ii) For every iy>l, the values in (4.2)~(4.4) are never attained, for any {X,Y) G F. 


Consider the corresponding solutions u'^{t,x) of (1.1), with graph 

{(u"(X,y), r(X,y), x'^(X,Y)); (X,y)GM 2 | ^ 

For t = 0, by (5.8) the corresponding sequence of initial values satisfies 

h^(0,-)-^ollcUMD = HiO,-) -ui\\cH[a,b]) = 

for every bounded interval [a, b]. 

Next, consider a cutoff function rj G such that 

r]{x) = 1 if |x| < r, 

r]{x) =0 if |x| > r + 1, 

with r » p sufficiently large. For every v >1, consider the initial data 

+ (1 - v)uo , u\ = rjui + (1 - r])ui. 


(5.9) 


(5.10) 


By (5.9) we have 


= 0. (5.11) 

Moreover, if r > 0 was chosen large enough, we have 

iY{t,x) = vy{t,x) for all {t,x) G [0,T] x [—p,p\, 

while iX remains on the outer domain (5.4). The above implies {uq,u'() G B for all > 1 
sufficiently large, proving that "D is dense on Bs- 


lim ||(u[(-uo, 
u—yoo 


5. To complete the proof we need to show that, for every initial data (uq, ui) G V, the solution 
u{t,x) of (1.1) is piecewise on the domain [0,r] x R. 

By the previous arguments, we already know that u is on the outer domain (5.4). It thus 
remains to study the singularities of u on the inner domain [0,T] x [—p, p]. For this purpose, 
call {u,w,z,p,q,t,x){X,Y) the corresponding solution of (2.16)-(2.20), with boundary data 
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as in (2.21). By (5.6), every point of the inner domain is contained in the image of the square 
r in (4.1). 

Consider a point (XojYo) £ h. Two cases can occur. 

CASE 1: w{Xo,Yq) ^ tt and z{Xq,Yq) ^ tt. By (2.19)-(2.20) it follows 

det r ^ _ (l + cosrc)p (l + cos 2 ;)p (l + cosz)^' (l + cosr(;)p ^ ^ 

\ tx ty J 4 4c 4 4c 

Hence the map {X,Y) i— {x,t) is locally invertible in a neighborhood of (Ao,lo). We can 
thus conclude that the function u is in a neighborhood of the point {t{Xo,YQ), x{Xq,Yq)). 

CASE 2: rc(Ao,yo) = In this case we have either wx{Xo,Yq) ^ 0, or else by (2.17) 

(y\ 

wy{Xo,Yo) = \ {cosz + l)q ^ 0. (5.12) 

OC^(U) 

Indeed, we always have c{u) > 0 and q > 0. Moreover, by construction the values {w, z, wx) = 
(tTjTt, 0) and {w,wx,c'{u)) = (7r,0,0) are never attained in E. This implies (5.12). 

By the implicit function theorem, we thus conclude that the sets 

= {(A,y)Gr; w{X,Y) =Tr}, = {(A,y)Gr; z{X,Y) = Tr} (5.13) 

are the union of finitely many curves. 

The set of points {t,x) where u is singular coincides with the image of the two sets S'^, 
under the map 

(A,y) ^ A(A,y) = {t{X,Y), x{X,Y)). 


6. To complete the proof, we study in more detail the images of the singular sets S'^, S^. 

By (4.5) there can be only finitely many points inside T where w = n and wx = 0, say 
Pi = {Xi,Yi), i = 1,... ,m. Moreover, by (4.6), at a point (Aq;I^ o) G 5'"' D 5^ we have 

wx / 0, Wy = 0, Zx = 0, zy 0. 

Therefore, as shown in Eig. 2, the two curves {tc = vr} and {z = vr} intersect perpendicularly. 
As a consequence, inside the compact set T, there can be only finitely many such intersection 
points, say Qi = (A-, A/), f = 1,..., n. 

After removing these finitely many points Pi,Qi, we can thus write 5“ as a finite union of 
curves 'jj of the form 

7,- = {(A,y); X = <Pj{Y), ai<Y<b,}. (5.14) 

for suitable functions 7 j of class C^. We claim that the image of A( 7 j) is a curve in the t-x 
plane. To prove this, it suffices to show that, on the open interval ]aj,bj[ the differential of 
the map 

y ^ (x(,^,-(y),y), t(</>,-(y),y)) 

does not vanish. This is true because, by (2.20) 

^t{cl)j{Y),Y) = tx-4>'j+ty = + > 0- 
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Indeed, 2: 7 ^ tt while c{u), q > 0. 


As shown in Fig. 3, restricted to the inner domain [0, T] x [—p, p] in the t-x plane, the singular 
set A{S^) is thus the union of the finitely many points 

Pi = A(-Pi), i = 

Qi = A{Qi), i = 1,... ,n. 


together with finitely many curves Fj = A(7j). The same representation is valid for the 
image A(5'^). This concludes the proof of Theorem 1. □ 


6 One-parameter families of solutions 


In this section we study families of conservative solutions u = u{t, x, A) of (1.1) depending on 
an additional parameter A € [0,1]. We thus consider a 1-parameter family of initial data 

u{0,x,X) = uo{x,X), ut{0,x,X) = ui{x,X), (6-1) 

smoothly depending on the additional parameter A G [0,1]. More precisely, these paths of 
initial data will lie in the space 

A = (^C3([0,l] X M) nL°°([0,l];i/^(M))) X (^C2([0,l] X M) nL°°([0,l]; l2(R))). (6.2) 

In particular, the map {x, A) i-7> uo{x, A) is three times continuously differentiable and the 
norm of no(-,A) is uniformly bounded for all A. Moreover, the map (x,A) i-)- mi(x,A) is two 
times continuously differentiable and the norm of tti('. A) is uniformly bounded for all A. 

By an adaptation of the previous arguments one obtains 


Theorem 2. Let the wave speed c{u) satisfy the assumptions (A) and let T > 0 be given. 
Then, for any 1-parameter family of initial data {uo,ui) G A and any e > 0, there exists a 
perturbed family (x,A) i-;- {uq,ui){x,X) such that 


{uo -uo, Ui-Ui 


< e, 


(6.3) 


and moreover the following holds. For all except at most finitely many X G [0,1], the conser¬ 
vative solution u = u{t,x;X) of (1.1) is smooth in the complement of finitely many points Pi 
and finitely many curves •jj in the domain [0, T] x M. 


Toward a proof, we shall need 


Lemma 6. Let the function u i-)- c{u) satisfy the assumptions (A), and let any M > 0 be 
given. 

Then there exists a dense set of paths of initial data T> C X such that, if{x, A) i-)- {uo,ui){x, A) 
lies in T>, then the corresponding solutions {t,x,u,w, z,p,q) of (2.17)-(2.20) with boundary 
data as in (2.21) have the following properties. On the domain F in (4-1) one has 
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{w,wx-,wxx) is transversal to the point (7r,0,0). 
(z, ZY, zyy) is transversal to the point (vr, 0,0). 
(w,z,wx) is transversal to the point (tTjTt, 0). 
{w,z,zy) is transversal to the point (vr,7r,0). 
{w,wx, c'{u)) is transversal to the point (vr, 0,0). 
(z, zy,c'{u)) is transversal to the point (vr, 0,0). 


Proof of Lemma 6. Consider any point (^O; Yo, Aq). Then, there exist 3-parameter families 
of perturbed initial data {uq,u^), 0 G such that the properties (1)—(3) in Lemma 4 hold. 
Indeed, it suffices to repeat all the arguments in the proof of Lemma 4 regarding Aq as a 
constant. For a fixed A = Aq, the perturbations in (3.12) are thus functions of s only, constant 

w.r.t. A. 

Combining these perturbations, as in the proof of Lemma 5, we obtain a map {X, Y, A, 9) i—>■ 
{u, w, z,p, q) for which all transversality conditions (i)-(vi) are satisfied. By the transversality 
theorem, for a.e. 6 the corresponding map {X,Y,X) i-)- {u^,w^ , z^,p^,q^) satisfies the same 
transversality conditions. This achieves the proof. □ 

Proof of Theorem 2. As in the proof of Theorem 1, we first choose p large enough so that 
all our solutions will be for {t,x) in the outer domain (5.4). 

For each A G [0,1], we denote by (u,w, z,p,q,x,t)(X,Y, X) the corresponding solution of the 
semilinear system (2.16)-(2.20). We choose M sufficiently large such that, for all A G [0,1], 
the inner domain [0, T] x [—p, p] is contained in the image 

A^(F) = {(t(A,y,A), x(A,y,A)); |A| + |y| < m}. 


By performing an arbitrarily small perturbation of the initial path of solutions we obtain a 
second path A u{-,X) such that, in the corresponding solution {u,w, z,p,q,x,t){X,Y, X), 
the transversality relations (i)-(vi) in Lemma 6 hold. 

Since the variables {X,Y,X) G F x [0,1] range on a compact, three dimensional set, this 
implies that the values in (i)-(vi) are attained only at finitely many points, say {Xi,Yi, Xi), 
i = 1,... ,n. Hence, for A ^ {Ai,..., A^}, the solution (t, x, u, w, z,p, q){-, •, A) does not attain 
any of the values in (i)-(vi), for {X,Y) G F. As shown in steps 5-6 in the proof of Theorem 
1, the corresponding solution u = u{t,x; A) is then piecewise smooth on the inner domain 
[0,T] X [-/),/>]. □ 

Remark 5. For a given solution u = u{t,x), define its singular set as 

= {(t,a:); u is not on any neighborhood of (t, x)}. 

In the above construction, one can regard Ai,..., An as bifureation values, where the structure 
of the singular set changes (Fig. 4). On the other hand, for A ^ {Ai,...,A„} the solution 
u{-,-;X) is strueturally stable. A small perturbation of the initial data does not change the 
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topology of the singular set. Based on the present analysis, we speculate that a theory of 
generic structural stability and a global classification of solutions to (1.1) can be developed, 
in analogy to the classical theory for ODEs [1, 16]. 





Figure 4: The singular set for a solution u{t,x;X). When the parameter A crosses one of the critical 
values Xi, the topology of the singular set changes. 
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